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Abstract 

In this paper we propose an extension of the Cahn method [1] to binary mixtures and 
study the problem of wetting near a two-phase critical point without any assumption 
on the form of intermolecular potentials. A comparison between Cahn's method and 
later works by Sullivan [2,3], Evans et al [4,5] is made. By using an expression of the 
energy of interaction between solid surface and liquids proposed recently by Gouin 
[6], we obtain the equations of density profiles and the boundary conditions on a 
solid surface. In the case of a convex free-energy, a one-dimensional solution of a 
linear problem is proposed for the density profiles between a bulk and on a solid 
wall. A non-linear model of binary mixtures [7] extending Cahn's results for simple 
fluids is also studied. For the case of a purely attractive wall we have established 
a criterion of a first order transition in terms of the structure of the level set of 
the homogeneous part of the free energy. Additively, explicit expressions of density 
profiles near the wall are proposed. They allow one to consider the adsorption of 
mixture components by a solid wall. 
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1 Introduction 



In 1977, Cahn [1] gave simple illuminating arguments to describe the inter- 
action between solids and liquids. His model was based on a generalized van 
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der Waals theory of fluids treated as attracting hard spheres [7]. It entailed 
assigning to the solid surface an energy that was a functional of the liquid 
density "at the surface". Three hypotheses are implicit in Cahn's picture for 
simple fluids: 

(i) In order for the liquid density to be a smooth function of the distance from 
the solid surface, that surface is assumed to be flat on the scale of molecular 
sizes and the correlation length is assumed to be greater than intermolecular 
distances (this is the case, for example, when the temperature T is not far 
from the critical temperature T c ). 

(ii) The forces between solid and liquid are of short range and can be described 
simply by adding a special energy at the solid surface. 

(Hi) The fluid is considered in the framework of a mean field theory. This 
means, in particular that the free energy of the fluid is a classical so-called 
"square-gradient functional". 

After Cahn, the problem of adsorption and wetting was studied by a statisti- 
cal method by Sullivan [2,3], Evans et al [4,5], respectively for gas and binary 
fluid mixtures. From the point of view of Sullivan and Evans et al one may 
view Cahn's approach as open to criticism for several reasons: 

(a) Cahn's treatment is based on phenomenological "square-gradient" version 
of van der Waals theory, which in contrast to the approach initiated by van 
Kampen [8] does not attempt to relate directly the properties of the non- 
uniform fluid to the interactions occurring on a molecular level. 

(b) The density adjacent to the wall vary strongly over the range of inter- 
molecular forces, consequently the gradient expansion approximation used in 
deriving the square-gradient theory is no longer valid. 

(c) Cahn's theory leaves unspecified a contribution due to the fluid-solid in- 
terfacial free energy. 

Evans et al [4,5] following Sullivan's approach [2,3] for simple fluids consider 
the special case of a contact between a two-component mixture near "the criti- 
cal end point" and a wall. They used Sullivan's grand potential to describe the 
solid-fluid and fluid-fluid interactions and tried to solve directly the problem of 
repartition of densities in a liquid (gas). Evans et al obtain a coupled system of 
integral equations for chemical potentials (cf. Eq. (6) in [4]). Then, to solve the 
system, it is necessary to know the interaction potentials between components 
and between solid wall and components: Evans et al assume an exponential 
interaction both for component-component and solid-components (as in [2,3]). 
Only such a hypothesis allows one to obtain two differential equations instead 
of the two integral equations (Eq. (10) in [4]). This assumption cannot be 
obviously valid for large classes of mixtures. Moreover, a special hypothesis 
(mixing rule) concerning interactions between components is assumed. Then, 
the mixing rule and exponential dependence allow one to obtain both the lin- 
ear relation between potentials and boundary conditions and the problem is 
reduced to the problem of the contact of one-component fluid with a wall. 
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The phenomenological "square-gradient" model is proposed in case of an in- 
finite non-homogeneous fluid or a fluid mixture as a small-gradient approx- 
imation by Widom [9] and Fleming et al [10]. The method is extended in 
mean-field theory for semi-infinite media in contact with a wall: as proved in 
[6], the fact that the densities are discontinuous at the solid wall does not 
disqualify the procedure used by Widom and Fleming et al and Cahn's treat- 
ment is valid for fluids and fluid mixtures near a critical point in contact with 
a wall. 

In this paper, we use the expression of a surface energy. The surface is assumed 
to be solid and interactions between solid and fluids are sufficiently short- 
range. The contribution of fluids is represented by a surface free energy with a 
density of the form Es(pis, P2s)> where p is and p 2 s are the limiting densities 
of the fluid components at the surface. The expression of the surface energy 
obtained in [6] is in the form: 

E S = '111 PlS ~ 721 P2S + \ (7l2 Pis + 722 Pis + 2 7 3 2 P1SP2S) ■ ( X ) 

This expression represents first terms of a more complex expansion. It is an 
extension with explicit calculations of the widely known expression due to 
Nakanishi and Fisher [11] and examined in a review paper by de Gennes [12]. 
All the coefficients 7^ can be calculated explicitly after the particular form of 
interaction potentials was chosen. For example, in the case of London forces, 
the values of coefficients related to the densities of the two fluids at the surface 
are [6] 

Pin 

7ll " ^ 5 2 P3 , 721 - 12 6 2 Ps , 
kin k 2 ix 

712 " ^ ^2 > 722 - ^ S 2 , W 
k 3 TT 

732 - 24 

where p 3 is the density of the solid, i G {1,2} are the coefficients associated 
with intermolecular potentials of interaction between the fluids and the solid 
wall, ki, i G {1,2,3} are intermolecular potentials of interaction between 
the molecules of fluid % and themselves or between the two fluids and 8j = 

— (<7j+r), % G {1,2} are the minimal distances between the solid and molecules 

of the two species of the mixture, where <7j, % G {1,2} is the diameter of 
molecule of fluid i and r for the solid. Expression (1) allows us to estimate the 
influence of a solid wall on each component of a fluid mixture. Depending on 
the values of coefficients 7^, one can estimate the magnitude of the attraction 
or repulsion effects due to the wall. 

As our approach is also based on a mean-field approximation, we assume 
that variations of densities near the wall take into account several molecular 
ranges. Hence, it is possible to present the total free energy of the system 
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"fluids - wall" as the sum of a bulk free energy and a surface energy which is 
an additional contribution arising from the non- uniformity of the fluid near 
the wall. By using an extended variational principle, we obtain two boundary 
conditions at the wall and two partial differential equations for the density 
profiles of the components between a solid wall and a bulk. The complete 
set of boundary conditions and equations for densities allow us to obtain the 
profiles of densities in the following physical situations. The first is the study 
of the linear problem associated with the equilibrium of a two-component one- 
phase mixture near a critical point with a solid wall. The second is the study 
of the non-linear problem of the contact between a two-component two-phase 
mixture near a critical point and a wall. We get a condition of wetting and a 
first order wetting transition in terms of the level set of the homogeneous part 
of the free energy. 

To clarify the presentation some calculations are situated in Appendices. In 
Appendix 1, we present general calculations by using an extended variational 
principle applied to multi-component mixtures. In Appendix 2 we give an 
analytical representation of the profiles of densities connecting bulk and solid 
wall for a general form of the free energy of a two-component mixture near any 
critical point. These representations may be used to investigate the adsorption 
of fluid components of a mixture by a solid wall. 



2 Equations of density profiles and boundary conditions: general 



The general form of the free energy per unit volume of the mixture is proposed 
in the form [7,10,13] 



where V notes the gradient operator in the physical space V. The associated 
total free energy is 



The wall boundary 5 of D is endowed with a surface energy per unit area. 
The surface is solid and sharp on an atomic scale and the interactions between 
surface and fluids are sufficiently short range; the general form of the surface 
free energy per unit area used is 



results 




(3) 




E s = E s (p is ,p 2 s)- 



(4) 



Consequently, the free energy of S is 
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Then, the grand potential of the system "fluid mixture - wall" is 

£ = JJl EdV+ IL EsdS - 

The condition of extremum of the energy £ based on hypotheses (3) and (4) 
yields (see for details Appendix 1): 

- Equations of two profiles of component densities: 

BE 

where — — — - is the vector whose components are the partial derivatives of E 
d(V Pa ) 

with respect to the components of V p a and div is the divergence operator. 

- Two boundary conditions at the solid wall: 

dE dE s 

n— — + — ^ = 0, a = 1,2, 6 
#(VpJ dp a 

where n is the external unit normal vector to T>. 

Equations of equilibrium (5) are the same as in [7] given for the one-dimensional 
case. Conditions (6) generalize those proposed in [1]. 



3 The dynamical system associated with one-dimensional density 
profiles 



In the simplest case, the surface energy per unit area is given by (1) where 
the coefficients 7^ are expressed by means of a mean-field approximation 
through the potentials of the intermolecular interactions (see for example (2)) 
and the free energy per unit volume is of the form 

E = U( Pl ,p 2 ) + l - (Ci(V Pl ) 2 + 2D V Pl Vp 2 + C 2 (Vp 2 ) 2 ) , (7) 

where U (p ± , p 2 ) is the homogeneous free energy per unit volume and C\,Ci,D 
are constants such that the corresponding quadratic form is positive definite 
(we denote the free energy by U corresponding in [7] to —W). 

Let us consider the case of a flat plate wall defined by equation z = (see 
figure 1), where z denotes the one-dimensional coordinate orthogonal to the 
wall. The equations of equilibrium (6) associated with (7) are 
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Fig. 1. One- dimensional contact of a fluid with a wall. 



Cl d lP± + D d2 P* 



D 



dz 2 
d 2 Pi 
dz 2 



+ C 2 



dz 2 
d 2 p 2 
dz 2 



dU 

dpi 
dU 

dp 2 



(Pi,P 2 ) + ei 
(Pi,P 2 ) + e 2 



(8) 



where e± and e 2 are two constants of integration. 

These equations are complemented by the boundary conditions (6) at z = 0. 
By using expression (1) of the surface energy, we get 



dp 1 dp 2 

dz dz 

D dpi h c dp 2 

dz dz 



-In + IxiPi + 732P2 
-7 2 i + 732P1 + 7 22 P 2 



(9) 



We have to add the condition in the bulk (at z = +00) 



Pi = Ploo) P2 = P2c 



(10) 



4 Linear wetting problem 



We consider the case of a one-phase mixture (liquid or gas) in contact with a 
solid wall. The densities of the two- components and the temperature are close 
to critical conditions. Moreover, we assume that density variations are small 
enough with respect to bulk densities, i.e. 



Pi 00 
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such that we can consider a linearized problem associated with equations (8). 
Let us denote 




q 





A 




B 



-7 + r Pc 



/ d 2 U p 2 u \ 
~dpj d Pi d P2 
d 2 U d 2 U 



The matrix B is calculated in the bulk (pioo,P2oo)- Taking into account the 
definitions (11), we get the linearized problem associated with equations (8)- 
(10) in the form: 

d 2 r 

A U=* (12) 

dr 



A— = q + Tr at z = 

dz 



(13) 

r = at z = +oo. (14) 

The stability of the thermodynamic state of the bulk requires that the sym- 
metric matrix B is also positive definite. 

Let Xi> be the eigenvalues and the eigenvectors of the equation 

(B - xlA) h t = 0. 

Since B and A are symmetric and positive definite, \ 2 are positive. We can 
always suppose thatLj 

hiAhj =0, i^j. 
The solution of (12) satisfying the condition (14) is in the form 



r = 6jhj exp(— x { z) where x% > • 



(15) 



i=i 



Substituting expression (15) into condition (13), we get a linear system of 
algebraic equations for the unknown coefficients hi 

2 

Y,bi{r + x i A)h i = -< l 



which defines a unique solution 6, if det(hiThj + XihiAhj) ^ 0. In particular, if 
T is negligible (we assume that the wall is purely attractive), we get q = — 7 



1 ViVW = V (g> (JVW) denotes the bilinear form of vectors V and W with respect 
to matrix N; the bilinear form is symmetric when matrix N is symmetric. 
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(see (11)) and 



7 ik 



In such a case, the solutions satisfy conditions (13) and the density profiles 
fulfil the solution of linearized problem. Equations (15) yield different forms of 
density profiles. Depending on wall conditions, we may obtain both monotonic 
and non-monotonic profiles. This is similar to results of [7] in the non-linear 
case without a solid wall. In figure 2, we represent the different density profiles 
for each component of the mixture. We note that only one extremum point 
may appear for each density profile. This result is different from the results of 
Evans et al where density profiles are essentially monotonic. 




Fig. 2. In (a) and (b) the profiles of densities for a component are monotonic. In (c) 
and (d) the profiles admit an extremum of density for a component. The extremum 
is not on the wall but at some distance into the bulk measured on a molecular scale. 
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5 Wetting problem near a critical point for a two-component mix- 
ture 



5.1 The dynamical system 



In a two-phase region near a critical point at a given temperature T, the 
expression of the free energy per unit volume U associated with a phase equi- 
librium is of the form [7] 



U = al{{blx^ + yf + {y + tf). 



(16) 



The parameter t is an independently varied field characterizing the "distance" 
from the critical point (p lc , p 2c ), a and b are functions of the temperature. 
The variables x and y are defined through the transformation 







'a b\ 


r = AR, R = 


, A = 








^ c d J 









('■) 


Pc = 




w 




\P2cj 



(17) 



The scalars a, b, c, d associated with the physical properties of the mixture near 
the critical point depend on the temperature T. The constants of integration 
e± — and e 2 = t are already incorporated in U. With Eq. (17), the system 
(8) can be rewritten in the form 



AMA 



d 2 R 

~dl? 



V R C/, 



(18) 



where A* denotes the transpose matrix and Vr£7 means the gradient with 
respect to R. 

Following Rowlinson and Widom [7] we denote by 



M = AMA = 



771 xx 771 X y 



(19) 



Obviously, if A is positive definite, M is also positive definite, i.e. m xx > 0, 
rrixx m yy — m l y > 0. The boundary conditions (9) at the wall are 



M— = g + GR, where g = A*(- 7 + Tp c ) and G = A*TA . (20) 

dz 
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In the following, we choose in Eq. (16) a = l/\/2 and bo — 1 (to do this, we 
have only to change the values of coefficients of the matrix A defined by (17)). 

Hence, U = \ {{x 2 + yf + (y + t) 2 ) . 



The system (18)-(19) yields 



d 2 x , d2 y o / 2 , . 
m xx — + m xy — = 2x{x + y) 

§f §f (21) 

d x d z y 2 

mxy d? + myy d^ = X +2y + t 
System (21) admits the first integral 



- rn xx (^) 2 + 771 (^)(^) + - m (^Lf - U(x,y) = 0. (22) 
2 xxK dz J xyK dz JK dz J 2 yyK dz J V ,y; V ; 

This integral is similar to the integral of energy for mechanical problems. 

Substitution of boundary conditions (20) into the relation (22) yields necessary 
conditions for x, y at the solid wall. For simplicity, we consider only the case 
of an attractive wall (G is then negligeable) . Conditions (20) yield 

M- = g. (23) 

In fact, it is natural to expect that the results we obtain in the case of an at- 
tractive wall are closely similar to the results associated with the most general 
case. Relations (22) and (23) yield 

U(x,y) = k 2 , 

where k 2 = gM _1 g. Then, discussion of the wetting of a fluid mixture with a 
solid wall arises naturally from the drawing of the level curves of U(x,y) as a 
function of the parameter t. 



5.2 Connection between the dynamical system and Young's conditions 



For a solid wall in contact with phases a and (3, the contact angle 9 is de- 
fined with the help of surface free energies a along the solid surface (Young's 
conditions) 

a aP cosQ = a aS - crp S ■ (24) 

where the different subscripts designate phases adjoining the surface or inter- 
face. No value of 8 satisfies Eq. (24) unless 

o- a p > | cr aS - o-p S | . (25) 
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If the inequality (25) is not satisfied, one of the fluid phases completely wets 
the solid and there is no contact between solid and other fluid phase. In fact, 
the forbidden surface is replaced by a layer of the wetting phase and the surface 
free energy becomes the sum of two surfaces' free energies of the layer 

VaS = &a/3 + &/3S- (26) 

Condition (26) corresponds to the perfect wetting with the solid. The surface 
energies can be calculated by the formulas 

/+oo r +00 r +00 

(K+U)dz, v«s = J o (K+U)dz, vps = J Q (K+U)dz, 



where K = - m xx (-^-) 2 + m xv {^-){-^-) + - m vv {-^-) 2 and integrals are taken 

2 xxK dz J xyK dz JK dz J 2 yyK dz J B 
on different paths connecting phase a and phase (3 or a phase and the wall [7]. 



Let us note (ds) 2 = \- m xx (dx) 2 + m xy dxdy + ]- m yy (dyy 



From the first integral (22), we get 

r ( x i3,Vfj) 1 r {x a ,y<x) 1 

(T a p= I (2U(x,y))*ds, a aS = (2U{x,y))*ds, 

J (Xa,Va) J (xMcVMa) 



a ps = / (2U(x,y))*ds. (27) 

J (XMp,VMp) 

The integrals (27) are calculated on the paths associated with system (21) and 
the boundary conditions on the wall 

x = x Ma , y = y Ma or x = x Mp , y = y Mf3 , 

and in the bulks 

x = x a , y = y a or x = x p , y = y p . 



5. 3 Discussion of the wetting 



For a solid wall the value of k 2 is given. Hence, the discussion depends on the 
relative value of parameter t. 

(a) t > and large enough. 

In this case we are far enough from the critical conditions. In figure 3a the 
phases are in points A(x a ,y a ) and B{xp,yp). One obtains easily that x a = 
~Vt, y a = —t, xp = y/t, yp = —t. The points M a and Mp belong to two 
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different connected components of the level set U = k 2 . In the vicinity of A 
(or B), the energy U is a convex function of x, y and as in Section 4, it is 
possible to find the profiles of densities connecting A and M a or B and Mp, 
respectively. The integrals (27) are positive and a Q)/3 is large with respect to 
o" Qj 5 and vp,s- Then the relations 

&a,s < 0a,0 + &p,s and a p,s < &a,p + cr a ,s (28) 
hold and we are in the case of partial wetting with 9^0. 

(b) t > and small enough. 

This case corresponds to phases close enough to the critical point (see figure 
3b). The level set U = k 2 consists only of one connected component containing 
the points M a and Mp. The phases are at the points A(x a ,y a ) and B(xp,yp). 
They are very close with respect to the distance to the level curve. The super- 
ficial tension a a> p is small with respect to the free energies a aj s and &p,s- The 
values of a a> s and o~p t s are in general different and one of the two relations 
(28) is not satisfied. We are in the case where one of the two phases wets 
completely the solid wall. No contact appears between the other phase and 
the solid. For exemple, if relation (26) is satisfied, the phase (3 wets completely 
the wall. 



(c) t < 

The mixture has only one phase at the point A(0, — 1/2), which is the only 
singular point of the system (21). The energy U attains a minimum at the 
point A (we note that for t > this point corresponds to a saddle point, 
which is not associated with a bulk phase). The free energy U of the mixture 
is convex at the vicinity of A (figures 3c and 3d). If t is small enough, the linear 
solution for the profiles of densities obtained in Section 4 can be used. When 
t is large enough, the solution for the profiles of densities can be calculated 
analogously as in the Appendix 2. 



5.4 Some remarks on the profiles of densities 



The system (21) yields 



d 2 R ( 2x{x 2 + y) 
x 2 + 2y + t 



M u = L, ' (29) 



and admits the first integral (22): 



f M f _ V +v) >-<y + t) > = . (30) 
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a. 2 9.4 id) -i -o.S 



0.5 1 



Fig. 3. Drawing of level curves for the free energy U for different values of t. Points 
A (and B) correspond to the bulks. The bold curves are paths connecting a phase and 
the solid wall or two phases. The interfacial tension is calculated along these paths, 
(a) is the case of partial wetting with a non-zero Young angle. The other figures are 
associated with different cases of total wetting of one phase: in (b) with two phases 
and in (c) and (d) with one phase. 



When the densities x, y are far from critical conditions, t is negligeable with 
respect to x and y and (30) reads 

£"£-<* > +«'> a -»' = °- < 31 > 

d 2 V 

Let us denote V(z) = RMR. Then, by using (29) and (31), we get -— = 

(Jj/C 

3a; 4 + 5x 2 y + 3y 2 . The right-hand side is a positive definite quadratic form, 
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which implies that 



d 2 V 

dz 2 



> 0. 



Hence, 



If 



dV 

dz 



V(z) > V(0) + V'(0)z. 
> , it follows from here that V — > oo as z — > oo. The level curves 



2=0 



of V are represented in figure 4a. Hence, x or y must be an increasing 
function of z near the solid wall. For example, let x be an increasing function 
of z when z is small enough. Due to the fact that x — > ± as z — > +oo 
and t is small with respect to x(0), the representation of function 
of z has the form shown in figure 4b. Hence, the function x(z) is non- 
monotonic. In this case we may have also non-monotonic profiles of densities 
unlike in the treatment of Evans et al. Then, construction of an analytical 
solution may be done according to the algorithm proposed in Appendix 2. 



v(x,y) - V 2 > V ( 




Vt 



(b) 



Fig. 4. (a) represents level sets of V. Since x or y must be an increasing function of 
z near the wall, due to the limit conditions at infinity, x or y is a non-monotonic 
function of z. (b) is a representation of such a function. 



6 Conclusion 



Near critical conditions, by using a variational approach, we have obtained for 
an isothermal binary mixture in contact with a solid wall equations of equi- 
librium and boundary conditions which generalize those obtained by Cahn. 
With limit conditions in the bulk, they form a closed boundary value prob- 
lem. When the free energy of the mixture is a quadratic form with respect to 
the densities of components and their gradients, we get explicit profiles of the 
densities in the one-dimensional case. 

In the case of a purely attractive wall we have also established a criterion 
of a first order transition, when a contact angle against a solid wall becomes 
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zero. This criterion is formulated in terms of the level set of the function 
U(x,y): U(x,y) = k 2 , where k 2 depends on the boundary conditions. If the 
level set is a connected set, two multi-component layers exist: one layer with 
ordinary adsorption and the second one in contact with the wetting layer. If 
the level set is disconnected we have partial wetting. We have also shown that 
the profiles of density are typically non-monotonic. This is in agreement with 
Rowlinson and Widom [7] where infinite two-phase two-component mixtures 
where considered. 
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Appendix 1. Calculus of variations for fluid mixtures 

We study a two-fluid equilibrium, but the method can be extended to any 
number of components. The position of a two-fluid mixture is associated with 
two applications 

x = 4> a (X Q ), a = 1,2, 
where X Q denote the Lagrangian coordinates belonging to a reference space 
T> a associated with the crth component and x denotes the Eulerian co- 
ordinates in the physical space V [13]. The virtual motions of particles are 
deduced from the relation 

x = $ Q (X a ,5 a ), $ Q (X Q ,0) =0 a (X a ). 

Here e a , a = 1,2 are small parameters defined in a neighbourhood of zero. 
Virtual displacement £ a are defined by [13,14] 

<9<E> 

C =— i(X a ,e a ), a =1,2. (Al) 

de a k«=o 

At the solid boundary, the virtual displacement £ Q is subject to the conditions 

nC = 0, a = 1,2, (A2) 
where n is the unit normal vector to the boundary. 
Eulerian variations of densities are defined by 

6aPa = l[T Pa ^ £a \e a =Q and S <*Pp = Q > P^ a > a, 13 = 1,2. (A3) 

The variations (A3) are related to the virtual displacements (Al) by the for- 
mulae [14,15] 

S a p a = -div(p a Co) , <*= 1, 2. (AA) 
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The variations of the volume free energy are 

SaSv = J J JjaEdV, 

where 

dE dE 
5 a E = -T—5 a p a + f Q S a Vp a with f a = 



d Pa ^ a a Q! rot a d{ypa) . 

Since 

S a Vp a = V{5 a p a ) 

we get 

5 a £v = J J J v {jf 5 «P« + div M«/0) dV, (A5) 
5E 

where the variational derivative - — is defined by 

5E dE , 

1 — = a dlv f «- 

op a op a 

From relations (A2), (A4) and (A5) we obtain 

KSv = I / X ( P " V (17) C " " div (c«-^C + f« div (p a ( a )Y\ rff 

= / / L "« v © c « ^ - / L a f - div ^ )dS ■ 

The variations of the surface free energy are 

yy 5 -yy, ^f div(p « u ^- 

The grand potential of the system is £ = £v + £s an d its cuth variation is 
given by the formula 

s - £ = IIL "* v © c « ^ - //, ( n f « + S) dlv (p "° dS ' 



Denoting c„ = nf a + ^ S , we obtain 

<9p a 



I J s C a dw(p a C a )dS = jj s [ca dw s (p a U + c a (n^-^nj j d5 , 

where div,s denotes the surface divergence. Denoting by V5 the tangential 
gradient to S, we get finally 

s ° £ = I 1 L "» v © ^ + / L - ("If") ) ds 
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Consequently, the equations of equilibrium are 

_ / 5E \ 



0, a = 1,2, (A7) 



or 



dE / dE 

hv 



div 0/vT =0 a = 1,2 



and the boundary conditions are 

c a = 0, and V 5 (c Q ) = 0, a = 1, 2. 

Due to the fact that V^Cq,) = is a direct consequence of relation c a = on 
the surface, the only effective boundary conditions are c a = , i.e. 

n— — — + -— = 0, a = l,2. 

9(Vp a ) d Pa 



Appendix 2. Analytical representation of the profiles of densities 
of a two-component mixture for the wetting problem near a critical 
point 

Our purpose is to express analytically the profiles of densities of a two-phase 
mixture in contact with a solid wall. The free energy is given by (16) and the 
dynamical system for the profiles are given by 



d 2 x 



d 2 y 



m xx — + m xy — = 2x(x +y) 

it k. 



d x 



d y 



m *y dz 2 + m yy dz 2 



x l + 2y + t 



(A8) 



We notice that the matrix 



d 2 U 
dR 2 



6x 2 + 2y 2x 
2x 



is positive definite in the bulk phases (±xqq, yoo) where yoo — ~ x oo an d t = 
Xq 7^ 0. As in the Section 4, we can determine the positive eigenvalues x 2 
defined from the equation 

det (m? ~ ** M ) - °- 

where all the matrix coefficients are calculated in the bulk. We obtain 



X [m xx m 



yy ' m xy. 



+ x 2 (^x om X y - 4:xl m yy 



2m xx ) + 4x 2 00 = 0. 
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For | xoo | small enough, we get the two eigenvalues in the form 



and 



2 2 2 

V = X 2 = N 



Am 



j^-Xqq + C(|a; o| 2 ), where $ = 2m xx /D 2 , 



D 2 = m xx m yy - m 2 xy . 



(A10) 



We are looking for the solution of the system (A8) which goes to the equilib- 
rium states (±Vi, —t), t > at z = +oo in the following form 



x = x km exp ( - (Xk + fim)z) 



k,m>0 



Y Vkm exp ( - (Xk + fim)zj 



(All) 



k,m>0 



We assume that this expansion is valid for all positive values of z. We will 
show that this solution represents a two-parameter family. The values of the 
parameters will come from the boundary conditions (20). Substituting rela- 
tions (All) into (A8) and denoting S = exp(— Xz) and e = exp(-fiz), we 
get 

(AA; + nm) 2 5 k e m (m xx x km + m xy y km ) 

k,m>0 

= 2 x km xipXq r 5 q E +p+ + 2 x km yip5 e +p 

k,m,l,p,q,r>0 k,m,l,p>0 

and 

Y (Xk + fim) 2 5 k e m (m xy x km + m yy y km ) 

k,m>0 

= £ x km x lp 5 k+l e m+p + 2 £ y km 5 k e m + t. 

k,m,l,p>0 k,m>0 

The identification of terms S j e k yields 



xqo 



-t 



l o . 



6 L e 



XlO 



\yio 



= ciohio with hio 



+ O(\x 00 \) , 



V 



<$V : | " I = c 01 h 01 with h i = I " | + O(\x 00 \) 



x i 



— m 



xy 
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The vectors h i and h 10 are eigenvectors corresponding to eigenvalues A 2 and 
fi 2 defined by (A9) and (A10), respectively. The constants c w and c 01 are 
multipliers to be defined. 

In the same way, the terms associated with 5 1 e 1 , 5°e 2 and 5 2 e° are 



coiciohn with hn 



m 



xy 



D\x 00 



■hoi + 0(1) , 




= with h. 



m 



( ~ (rn xx 



02 



xy 



m yy + m l y ) \ 



2m, 



m 



xy 



+ O(\x 00 \) , 




= c 2 h. 20 with h 20 = 




The expansion of R truncated to the second order with respect to c w and c i 
is 

R = R 00 + ciohio exp(-Az) + c ih i exp(-(xz) + c iCi hn exp(-(A + fx)z) 

+00x1102 exp(—2fxz) + c 2 h 2 o exp(-2A^) + ... (-412) 

Because the solution R must be bounded as |xoo| goes to zero for all positive 
values of z, cio is at least of order |xoo| - Then, we may introduce a constant 
cio by the formula 



cio = .DIxqoIcio, D 2 = m xx m. 



yy ' m xy • 



Since the terms R o and c 10 h 10 are of the order of C(|a;oo|), we get in the 
vicinity of z = 



R = a ih i exp(-fxz) + c 2 10 h 2 exp(-2fxz) + O(\x 00 \), 



(A13) 



where all the coefficients are finite as |xoo| goes to zero and a i = c oi — 
mxyCoiCw- In fact, the term R o is negligible at the vicinity of the wall, but 
not in the bulk. 



For the sake of simplicity, we exhibit the boundary conditions only in the case 
of a purely attractive wall (when T is negligible). 
Then, the condition (20) on the wall reads 

M— = -A*7 at z = 0. 

dz 
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Or, by using (A13), we get 



/i M(a 01 h 01 + 2c 2 10 h 02 ) = A* 7 . (AU) 

By multiplying (A14) by h 10 and h 01 , we get a system of two scalar equations. 
By taking into account the equality h 10 Mh 01 = 0, the vector equation (A14) 
is equivalent to the system of two scalar equations 

2 / u c 2 h 10 Mh 02 = h 10 A* 7 (A15) 

/i (a 01 h 01 Mh i + 2c 2 h 01 Mh 02 ) = h 01 A* 7 . (A1Q) 

Eq (A15) defines c 2 and then, Eq. (A16) defines aoi- Hence the solution 
near the wall (A13) is completely determined. Consequently, we are able to 
determine the effect of the solid wall on the adsorption of each component. 
Due to the fact that the expansion (A13) depends only on Coi 2 and not coi, in 
this approximation we do not need to satisfy the inequality c i 2 > 0. 
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